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Abstract 

We present a new mechanism for B — > rfK , and analyze it strictly within perturba- 
tive QCD. We find that such a new mechanism may dominate non-leptonic -B-decays 
to light mesons. Within a reasonable parameter space, our prediction is in good agree- 
ment with the recent CLEO data on B^ — > rj'K^. We conclude that there is no room 
for an abnormally large b — > sg vertex from physics beyond the Standard Model. 
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Recently, the CLEO collaboration has reported a large 77' yield in charmless B decays as 
follows JLJ, 

BTZ(B ± -> r)'X s ) = (6.2 ± 1.6 ± 1.3) x 10" 4 (2.0 < P v > < 2.7GeV) , (1) 

and a corresponding large exclusive branching fraction 

BTl{B ± -> t/K*) = (7.1^-s ± 0.9) x 1(T 5 . (2) 

It is another great experimental achievement in rare B decays since the measurements of 
B — > K*"Y and B — > X s 7 which involve the so-called QCD- and electroweak-penguins. Since 
then, theoretical investigations on these decays have appeared, offering several interesting 
interpretations of the data, both within and beyond the Standard Model (SM) [|3|. In 
particular, the decay mode (2) makes itself conspicuous due to its surprisingly large branching 
ratio. 

Now let us examine the theoretical status of the estimate of the exclusive branching fraction 
B{B ± — > rj'K^). The standard theoretical framework to study non-leptonic B decays is based 
on the effective Hamiltonian which describes the decays at quark level, 

G 

H eff = — ^ X] V CKM C iOi , (3) 

and the BSW model || to estimate the hadronic matrix element < M\M 2 \ Qi\B >. 

An important feature of the BSW model in nonleptonic decays is the use of the factorization 
and spectator ansatz. It works remarkably well in the so called heavy-to-heavy transitions 
because the ansatz is consistent with the HQET, and therefore is almost certainly true. 
However, we think that it has never been tested, and so is much less justified, in heavy-to- 
light decays. 

Even given the validity of the BSW ansatz, yet another problem exits in estimating the 
hadronic matrix elements like 

< 7/ I u(l - 75)6 I B~ >, < K~ I s(l + 7s)« I >, 

< V I si5 s I >, < rf I uj 5 u I > . (4) 
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The Dirac equations of motion 



qq -- 

m q — Tflqi 

are used to get the elements for - 



-id^(q-f^ 5 q') 



(5) 



< T]' | M75M I > 

< K~ I s(l + 7 5 )m I > 
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m s +m u 
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With these relations, the amplitude for B ± — > rj'K ± reads 



Ml-m% F^\m\) f K 
a2 + aX Ml-m\Fr K -{m%)U 



+ (a 4 + 



2a 3 - 2a 5 + (a 3 - a 5 
2a 6 m|- 
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m s (mb — m s ) f^ w > 
M% -m^, F^\ml) f K 



(m s + m u )(m b - m u ) M 2 B - m\ F B ^ K ~ (to*,) 
< K~ I s Tm (1 - 75 )6 I fi- >< 77' I «7 M (1 - 75 )u | > 



} 



(9) 



Assuming the Dirac equation is valid for bounded fermions, the factor in the amplitude 
enhances the contribution of operator (9 6 , which is very difficult to understand. 

In an alternative way, we turn to estimate the amplitude using the perturbative QCD 
hard scattering formalism ||, which has also been extensively applied to 5-decays 0. The 
diagrams to be calculated are depicted in Fig. 1. However, we find that these amplitudes 
are rather small due to cancellations between the diagrams and the small b — > gg vertex 
@. Nevertheless, before one goes beyond the SM, the contributions in the SM should be 
examined carefully and exhausted. In what follows, we discuss a new type of mechanism for 
5-decays to light mesons in detail. 

The new mechanism is depicted in Fig. 2. This mechanism is motivated by the fact that 
both the recoil between 7/ and K~ and the energy released in the process are large. The 
gluon from b — > sg vertex carries energy about Ms/2 and then materializes to rj' and emits 
another hard gluon to balance color and momentum. The momentum squares of the gluons 



scale as ocM B 

k\ = ((1 - x)P B - yP K f Ml - (1 - x)y(M 2 - M 2 ,) , 
and*; 2 = (xP B - (1 - y)P^) 2 = -x(l - y)M% , (10) 

where x and ?/ are the momentum fractions carried by the collinear quarks shown in Fig. 2. 
For self-consistency of the co-linear picture used here, the terms ~ x 2 M B are neglected since 
they are at the same level of the transverse momentum square of the quarks in the bound 
state B meson. Using mean values < y >~ |, < x >~ e B with e B ~ 0.05 — 0.1 ||, we get 
< k\ >~ 12 GeV 2 and < k\ >~ 1 GeV 2 , which are large enough to justify a perturbative 
calculation. 

The soft contributions are parameterized in terms of wave functions of the bound states. 
In the spirit of Ref . || , we neglect the transverse components of quarks and take the wave 
functions for B~ and K~ as 

^ = 7i %MvH 4k , (ii) 

where 1q is an identity in color space. In QCD, the integration of the distribution amplitude 
is related to the meson decay constant 

J <f) K {y)dy = -^mfK , J 4> B (x)dx = ^=/b • (12) 
We write down the amplitude of Fig. 2 as 



KA f, A , i ^ Tr bs qf r M ( j/+ M B ) lblu ] Ae^k la k 20 C eff C F g 3 s 

M = J dxdyM*) 7m ^ i - h)klkl My) 



(13) 



where r M is the effective b — > sg vertex known for years 
G F 9s 



K = %7lW{^7,. - k!i)L - F^a, v k\m b R} . (14) 



We have used a Lorentz and gauge invariant amplitude < g a (ki,€i)g b (k2,€2) \ rf > in Eq. 
(0), with the shorthand notation G^ v ~ k^e u — e M fc„, and it reads P,[10[| 

< <?V | rf >= 5 ab A v ,F(kl k\)^°G\ v G 2 pa (15) 



with A v > = 2C % ff and the form factor F(kf, = 2k v ' k , which will be extracted form the 
J/ty decays. Cf is a color factor Cf = -^^Tr[T a Tb]5 a b = |. 
Finally we obtain 

M = ^a 2 s C eff C F 32(V is V*) 



J dxdycf) 



,F{k 2 \p ■k 1 q-k 2 -p- k 2 q ■ fei] + F l 2 M B m h \q ■ k 2 k\ - q ■ k x k x ■ k 2 ] 
B ' X)MV) k^kM > ( 16 ) 



where the momenta k%, k 2 , p, q can be read off from Fig. 2. We note that from Eqs. (10,16) 
the singularities of the gluons-propagators are located at the point y = 1, which is just the 
end point for the wave function ^k{u)- However, it is well known that the value of bound 
state wave function at its end point is exactly zero in QCD. 
We can extract C e ff form the data J/^f — > 7/7 to remove ambiguities: 
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r(J/fr _> rfj) 4«K)^ R 2 x- I H^jx) 
r(J/*^e+e-) ~ vr a e M|^ Vd ° e//j 54 



(17) 



where the x- \ H PS (x) | 2 can be found in 0. For x = 1 — -g^-, x- | H PS (x) | 2 = 54.8. We 
extract C e // — 0.075 GeV. In order to get quantitative estimates, we take the wave function 
as [| 

<t>B{x) = —^8{x - e B ) , and (j) K {y) = V3f K y(l - y) , (18) 

and t b = 1.62 ps, f B = 140 MeV, / K = 113 MeV (corresponding f n = 92 MeV) and 

I VtsVtb |= 0.044. For the strong coupling at the energy scale < k\ >, we respect the choice 

a s = 0.38 in Ref. 0. The main uncertainty of our estimate exists in the distribution function 

4> B (x). At present, it can not be derived from the first principle of QCD. However, it is easy 

to understand that the distribution function is peaked sharply at one point due to the heavy 

b quark carrying most of the momentum of the B meson. The peaking position is expected 
at x = Mm^M, ^ Aqco^ 

M b M b 

From Eq. 10, 16, 18, it is easy to see that the integration in Eq.lQ is infrared safe because 
the pole in ljk\ is canceled by the K meson wavefunction and k\ is always positive. To 
confirm that | M. \ is dominated by hard gluon contributions, we define R as the ratio of 



I M. | with cut on k\ to that without cut. Taking Aqqd — 200MeV, we list the numerical 
results of R in Table I to show that the integration in Eq.lQ is indeed dominated by hard 
gluon contributions. For somewhat complete investigation on this mechanism, the next to 
leading order contributions(high twist) should be included, however, it is complicated. We 
shall discuss it elsewhere. 

The numerical results of B1Z(B ± — > r/K^) are displayed in Fig. 3 as a function of the 
peaking position parameter e B . 

We can see that our predictions are in a good agreement with experimental results in the 
region of e B = 0.05 ~ 0.07. Furthermore, if the contributions of Eq. (9) estimated from the 
conventional way (which may contribute up to BTZ(B ± -> rfK*) = (1 ~ 4) x 10~ 5 P,P|) 
are taken into account, the SM predictions turn out to be in agreement with the CLEO data 
within the la level in the whole parameter range of e B = 0.05 ~ 0.1. We conclude that 
B1Z(B ± — > rj'K ) is not "surprisingly large", and the mechanism in the Standard Model 
presented here seems sufficient for explaining it. 

Note added: After finishing this paper, we became aware that a similar idea appeared in 
hep-ph/9710509| [jl3|| . But our physical picture, calculation method and conclusion are all 



different from theirs, which used a gluon- diffusion picture and calculated it by using effective 
Hamiltonian method. However, we here present a hard scattering picture and solve it on 
perturbative QCD. Acknowledgment We thank H.Y. Cheng and M. Drees for careful 
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Figure 1: Tzme and space H^e strong penguin diagrams for B ± — > r]'K ± in usual way. Dia- 
grams suppressed by V u b are neglected. Blobs depict K and rj' 
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Fig. 2 

Figure 2: New diagram for B ± — > r]'K ± . x, 1 — x, y and 1 — y are momentum fractions 
carried by the quarks. 
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Figure 3: B1Z(B — > r]'K ± ) as a fuction of e B . The horizontal thick solid lines show the 
CLEO measurement (with ±la error bar) and the thin curve is the contribution from the 
new mechanism presented here. 
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Cut 


R(Jb=.05) 


R(£b=.06) 


R(£b=.07) 




96% 


97% 


98% 


- k 2 2>4N<2CD 


76% 


87% 


95% 


- kj>9M>cD 


50% 


58% 


64% 



Table 1: Numerical results for R for different cuts and 6b ■ 
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